Full Counting Statistics of a charge shuttle 
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We study the charge transfer in a small grain oscillating between two leads. Coulomb blockade 
restricts the charge fluctuations in such a way that only zero or one additional electrons can sit on 
the grain. The system thus acts as a charge shuttle. We obtain the full counting statistics of charge 
transfer and discuss its behavior. For large oscillation amplitude the probability of transferring h 
electrons per cycle is strongly peaked around one. The peak is asymmetric since its form is controlled 
by different parameters for h > 1 and n < 1. Under certain conditions the systems behaves as if 
the effective charge is 1/2 of the elementary one. Knowledge of the counting statistics gives a new 
insight on the mechanism of charge transfer. 

PACS numbers: 73.23.Hk, 72.70.+m, 85.85.+J 



I. INTRODUCTION 

A few years ago Gorelik et al. showed that a small 
conducting grain can undergo a mechanical instability if 
it is trapped in a harmonic potential between two leads 
kept at a constant voltage bias (see Fig. QJA When 
the central grain is charged, the electrostatic force in- 
duced by the leads pushes it towards one of the elec- 
trodes, increasing the probability that the excess charge 
be discharged. Since the resistance depends exponen- 
tially on the distance, even small oscillations can largely 
amplify the probability of transmission. Excitation of 
this nanomechanical system at one of its resonating fre- 
quencies can be generated by the stochastic tunnelling of 
electrons from the leads. Since the charge state of the 
grain is correlated with its position, under certain condi- 
tions, the energy accumulated in the mechanical systems 
increases indefinitely, leading to an instability. The en- 
ergy pumped depends on the oscillation amplitude up to 
a maximum value determined by the number of charges 
that can accumulate in the grain at each cycle. After 
that point there is no additional gain in increasing the 
amplitude and the grain stabilizes at a fixed oscillation 
amplitude for which the energy pumped exactly balances 
the energy dissipated. This scenario has been investi- 
gated both in the incoheren1jii2i£ and in the quantum 
cage| 4 i 5 4 6 4 7 i 8 4 9 There are strong indications that Parks 
et al. have observed this phenomenon in Cgo molecules 
oscillating between two leads^ 

An other possibility to drive the oscillations is to use 
an external alternate electric field acting on a a can- 
tilever. The amplitude can thus be tuned independently 
of the source/drain voltage bias, at least in principle. 
This case has been experimentally realized by Erbe et 
al. who observed a current of 0.11 electrons per cycle at 
low temperature induced by the oscillation of the central 
grain. When the leads and the grain are superconduct- 
ing the existence of phase coherent transport as been 
proposedii2ii&ii 

Many papers studied theoretically the conditions for 
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FIG. 1: Simple schematization of the system. A small grain 
oscillating in an harmonic potential between two leads kept 
at a constant voltage bias. 



the realization of the instability or considered the de- 
pendence of the current on the external parameters. 
Only few papers investigated instead current fluctua- 
tions. Weiss and Zwcrgeui^ calculated the average num- 
ber of electrons transmitted and its fluctuation during 
a single cycle of a shuttle driven at a given frequency 
and amplitude. This quantity differs from the noise 
actually measured since typical measurement times are 
much longer than one period of oscillation. Correlations 
of charge fluctuations on different cycles are then im- 
portant, as we discuss in the following. Other authors 
considered the finite frequency noise in superconducting 
shuttles^ or the telegraph noise induced by the switch- 
ing between two mechanical modes in a two-oscillating- 
grains device^ A related problem is the fluctuation of 
the acoustoelectric current carried by surface acoustic 
waves propagating along a ballistic quantum channel*^ 
The full counting statistics (FCS) of charge transfer 
in a shuttle has not been considered so far. Recently, 
powerful techniques have been developed to calculate the 
probability than n electrons are transferred during a mea- 
surement time t in electronic devices jil^LSSiSI The FCS 
contains much more information on the dynamics of the 
charge transfer than the current or the noise alone. This 
will be particularly clear in this problem since few elec- 
trons are involved in the tunnelling, and the probabil- 
ity that h electrons per cycle are transmitted is actually 
a fundamental quantity. For a well developed shuttling 
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regime the noise to current ratio is expected to be small, 
since the number of electrons shuttled at every cycle is 
determined by the Coulomb blockade conditions and thus 
it does not fluctuate as it happens in a purely stochastic 
tunnelling. This implies that the probability distribu- 
tion has a small width, but its actual shape still depends 
on the physical parameters of the junction, like voltage 
bias, tunnelling probability, or oscillation amplitude of 
the shuttle. The importance of studying theoretically 
the FCS is thus twofold: first it is, at least in principle, a 
measurable quantity and secondly, its knowledge allows 
to infer detailed information on the mechanism of charge 
transfer. 

The paper is organized as follows. In Sec. [H]the tech- 
nique for the calculation is developed and the equations 
for the numerical approach are obtained. In Sec. IHII the 
FCS is calculated analytically for small and large oscil- 
lations. In Sec. I1VI the general (numerical) results are 
discussed and compared with the analytical ones. Sec- 
tion |V| gives our conclusions. 



II. FULL COUNTING STATISTICS FOR AN 
OSCILLATING GRAIN 

Our aim is to calculate the FCS of charge transfer in 
an oscillating grain between two leads. We assume that 
the charge transfer is described by the "orthodox theory" 
of Coulomb blockadei^ 2 * 2 ^ In this regime the tunnelling is 
incoherent and the dynamics is governed by a standard 
master equation. Within these assumptions Bagrets and 
Nazarov have developed an elegant and efficient tech- 
nique to derive the FCS for the static case^ We will use 
their method generalized to a moving grain. 

Since we arc interested to a single grain structure, the 
state of the system is completely determined by the prob- 
ability pk of having k additional electrons in the island. 
For simplicity we consider the case where the voltage bi- 
ases guarantee that only the two states, k = 0,1, are 
available, and that only two events are possible: either 
one electron jumps on the island from the left lead [with 
transition rate T^(t)] or one electron on the island (if 
present) jumps to the right lead [with transition rate 
r#(i)]. Within these assumptions the time evolution of 
the probability is given by a master equation. We write 
it following the notation of Ref. [2(1 
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\p(t)) = -L(t)\p(t)) 



(1) 



where the (classical) probability is represented by a state 
in a vector space: \p) = {po,Pi}, and L is the matrix 



L(f) 



r L (t) -r fl (t) 
-r L (t) r R (t) 



(2) 



In Ref. [2(J it is shown how the FCS can be obtained 
by calculating the time evolution of the probability with 
a modified operator L(t). The central quantity is P to (n), 



the probability that n electrons have been transmitted 
during a measurement time t a . This quantity is indepen- 
dent from the initial condition in the limit of large t D . 
From the technical point of view it is easier to calculate 
the generating function St„(x) : 



,s to (x) 



n=0 



(3) 



From S one can easily obtain all cumulants: n = 

d(-s)/d(i X )\ x =o, (n-ny = a 2 (-5)/a( lx ) 2 | x=0 , etc. 

Let us count electrons crossing, for instance, the left junc- 
tion. According to the prescriptions of Ref. |2{J the modi- 
fied operator L x (t) is obtained from Eq. J5J by multiply- 
ing the lower off-diagonal matrix element by the factor 
e lx . This factor keeps track of the electrons that cross 
the left junction during the time evolution. The gener- 
ating function St„ (x) IS then simply given by the formal 
integration of the modified time-evolution equation: 



<<z|Texp - / L x (t')dt' |p(0)>, (4) 



S to (x) 



where \p(0)) is the probability at time t — 0, \q) = {1, 1}, 
and Tcxp is the time ordered exponential. The derivation 
of this equality was done in Ref. |2?3 for the static case 
where T's do not depend on time. Following the steps of 
their proof it is not difficult to verify that Eq. holds 
also in the dynamical case of interest here. The main dif- 
ference is that for the static case (and the zero frequency 
noise) one can restrict to the study of the eigenvalues 
of L, since the time ordering becomes immaterial when 
L does not depend on time. In our case instead time- 
ordered exponential must be evaluated explicitly. 



A. Specific expressions for the shuttle 

Let us now consider explicitly the time dependence of 
the tunnelling rate. With good accuracy one can assume 
that the dependence of Tl /R on the position of the grain 
is exponential 



L/R 



r° exp{=Fa;/A} 



(5) 



(we assume = Tn = T° for x = 0). Here x is the 
shift of the grain from the equilibrium position and A is 
the tunnelling length (see also Fig. [JJ. We will consider 
the case of sinusoidal oscillations of the grain. This can 
be driven by an external device like in the experiment of 
Ref. 0, or it can be induced by the voltage bias between 
the left and right leads>i*i& In both cases 



IW*) =r°exp{Tasin(wt)}, 



(6) 



where a = x max /X is the dimcnsionlcss ratio of the os- 
cillation amplitude to the tunnelling length and u> is the 
frequency of oscillation. It is also convenient to rescale 
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the time by a; -1 and define <f> = uit. With this sub- 
stitution the problem is fully characterized by the two 
dimensionless parameters: a and T = T°/u>. From the 
physical point of view, T gives the probability that an 
electron in the static junction with x = can tunnel on 
or off the grain in the time 1/oj. We will see that the de- 
pendence on a of the FCS will be qualitatively different 
if r is smaller or larger than one. 

The interesting physical quantity is the FCS for a long 
measurement time t . We choose t Q to be a multiple of 
the period: t Q = 2nN/uj, with N integer. The FCS of 
charge transfer during N periods is then given by 

e -s N (x) = ( q \A N \p{0)) 



where 



and 



A = Texp j- J L x (</>')#' 



L x {4>) 



Pg-asin t 
Pg-a sin (j> 



_p g asin < 
p e asin<^ 



(7) 



(8) 



(9) 



In the case of interest of large N the FCS is given by the 
eigenvalue Xm{x) °f ^ that has the maximum absolute 
value: Sn(x) — ~ NIii[\m(x)]- From Sjv(x) one can 
calculate directly the particle current and noise reduced 
to a period: I = h = n/N and P = 2(n — ri) 2 /N. 

To obtain the probability of having transferred n elec- 
trons during N periods it suffices to invert Eq. 0: 



e -S N (x)-ixn _ 



For large N the saddle point approximation gives a very 
accurate estimate of this integral: 

]n[P N (n)}/N = ln[X M (xo)] - W (11) 

where n = n/N is the number of electrons transferred 
per cycle and xo satisfies the equation: 

1 dX M 



Am(xo) dx 



(12) 



X = X() 



We find that Eq. II 1211 is solved by xo pure imaginary. 

The problem is now reduced to the evaluation of the 
time-ordered product that enters the definition of A. 
This can be done numerically by integrating the system 
of differential equations 

^b(0> = -4W)li#)> (is) 

with the two initial conditions \p^'(cl) = 0)} = {1,0} and 
\p^ 2S> {4> = 0)) = {0,1}. One can readily verify that the 
matrix with columns \p^ x \(j) = 27r)) and \p^{4> = 2tt)) 
coincides with A. In the case of interest the numerical 
task is not hard, nevertheless discussion of tractable an- 
alytical limits greatly enhances the understanding of the 
results. We thus discuss in the next section the small 
and large amplitude limits before presenting the numer- 
ical results for the general case in Sec. IIVI 



III. ANALYTICAL LIMITS FOR SMALL AND 
LARGE AMPLITUDE 

A. Static case and fractional charge 

For a = we have a standard static single electron 
transistor. Since L x does not depend on <j>, the time- 
ordered exponential becomes a simple exponential 



-2ttL v 



(14) 



and the generating function for large N can be obtained 
by diagonalization of the L x matrix. The smallest eigen- 
value in modulus gives £W(x) aTi the leading order. The 
generating function reads 

-5 A r( X )/^-27rr(e' i x/ 2 -l) , (15) 

in agreement with the result obtained with a different 
technique in Ref . I24L The current and the noise are thus: 
I = P = nT with a Fano factor F = P/2I equal to 1/2. 

Even if Eq. I|15f) has been derived before its meaning 
has not been fully discussed and it is worth a short digres- 
sion. Tunnelling through a single barrier is a Poissonian 
process. The generating function and the probability in 
this case is 



S to (x)=n(e tx -l) and P t » = e""-^ , (16) 



with fi the average number of charges transmitted during 
the time t Q . A general feature of the generating function 
is the 2-7T periodicity. It is a manifestation of the discrete 
nature of the charge and follows directly from the defini- 
tion 121 . It is thus surprising that the generating function 
(|15fl is periodic of An as if the elementary charge was not 
one, but 1/2. This happens only when Tl = and 
for long measurement times (N — > oo), in all other cases 
S(x) is periodic of 2ir. For instance if Tl ^ Tn, even 
in the large N limit we have the following 27r-periodic 
generating function 



Sn(x) 
N{2tt/lo) 



(17) 

The result (|15J) for Tl = Tr indicates that the charge 
transfer in our system (Fig. [21 a) is equivalent to that 
happening in a tunnel junction with charges 1/2 emitted 
with probability T (Fig. b). 

This can be understood with simple arguments. Let 
us consider a sequence of events in our system. These 
can be of two types, either tunnelling from the left lead 
to the grain (type L), or tunnelling from the grain to 
the right lead (type R). Since the rates for both events 
are the same (r°), the system has the same probability 
per unit time to switch to the other state. The statistics 
of the switching events (i.e. that either R or L occurs, 
without specifying which one) follows thus a Poissonian 
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(a) 



result is 




(b) 



e/2 



FIG. 2: The two equivalent systems for the statistics of charge 
transferred in the long measurement time limit: (a) single 
electrons transistor with equal probability F of hopping from 
the left lead to the grain and from the grain to the right 
lead (b) a fictitious system of charges 1/2 that tunnel with 
probability T. 



distribution, since all events are independent (rates do 
not depend on the initial states). 

To obtain the statistics of charge transfer from the 
statistics of switching events it is enough to remember 
that every two switching events one charge is transmit- 
ted. We can thus associate the transmission of a fictitious 
1/2 charge at each switching event. It is clear that the 
counting statistics of the fictitious charge coincides with 
that of the true charge, apart from a possible charge 1/2 
mismatch that is irrelevant for long measurement times 
(t o T — > oo). Thus we proved that the counting statistics 
of our system coincides with that predicted for a tun- 
nelling junction of charges 1/2 transmitted with rate T°: 
simple classical correlations can induce current fluctua- 
tions typical for fractional charges! A similar "fractional" 
behavior was also found and discussed by Andreev and 
Mishchcnko for charge pump in the Coulomb blockade 



regime 
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One should keep in mind that for any finite measure- 
ment time the periodicity of the generating function re- 
mains 27T, it is only the leading term for t a — > oo that 
is Air periodic. We thus expect that the Fourier series 
is not uniformly convergent and that for any finite t Q 
higher moments will definitely depart from the predic- 
tion obtained with Eq. (|15|) . This can be verified by 
calculating the first correction to Eq. I) loll (we recall that 
2tt NT = t T ): 



S 



N 



2ttAT 



= e^ 2 - 1 



1 



27TAT 



log(l + cosx/2) + 



that holds with accuracy 0(e' 



-2irNF 



(18) 

). From Eq. JTSJ we 



derive all moments M q = d q {—S)/d{ix) q - For q even the 



M„ 



2-kNT 
2i 



(-1)9/2+1 4C(g) 

2ttAT tt? 



(l-2-«)(g-l)! 



(19) 

while for q odd the correction in Eq. (|18fl gives a vanishing 
contribution. For large q and fixed N the second term 
in Eq. {HI is approximatively 4(g/e7r) 9 /(2nNT). This 
means that for q large enough M q will depart from the 
prediction of Eq. I|15|) . This is the way in which the sys- 
tem may reveal the true nature of the elementary charge. 
Either by a short time measurement of the first moments, 
or by the long time measurement of higher moments. 

To explain while for asymmetric tunnelling rates the 
periodicity is 2n even for large t D it is enough to no- 
tice that the above discussed mapping cannot be real- 
ized when ^ Tr. It is worth mentioning that for not 
too large asymmetries the fractional charge remains mea- 
surable, as it is clear from the dependence of the Fano 
factor^ 



F 



r 2 « 



(20) 



But again one expects that departure from the prediction 
of Eq. (|15f) will increase with the order q of the moment, 
and for any small asymmetry it will become large for 
q large enough. Generating functions with periodicity 
induced by smaller fractions of the elementary charge 
can be obtained with several islands. 



B. Large oscillation amplitude 

Let us now consider the opposite limit of large oscil- 
lation amplitude of the shuttle. In this limit, since for 
most of the time the ratio Tl/Tr is either very large or 
very small, we can assume that (i) for < 4> < ir the 
quantity Tl vanishes identically and (ii) for n < (f> < 2ix 
the opposite holds: Tr = 0. The approximation becomes 
exact for r C 1, since in that case electrons can tunnel 
only when the shuttle is near one of the two leads. 

Within this approximation A can be obtained analyt- 
ically. As a matter of facts in region (i) po{4>) + Pi(.4>) is 
conserved, since the matrix element L that multiplies e %x 
vanishes. (We recall that pk(4>) is the probability that k 
electrons are present in the grain at time (p.) The intro- 
duction of the counting field normally breaks the conser- 
vation of the probability Po{4>) + Pi{<P) = const. Using 
this conservation and integrating the remaining differen- 
tial equation, we obtain for region (z) 



Po(tt) 
Pi(tt) 



Po(0) + (1 



a)pi(0) 



(21) 



where 1 — a is the probability of transferring one electron 
during half cycle given with 



exp 



3 asin<* 



(22) 



5 



In region (ii) po(<f>) + e lx p\(<j)) is conserved and wc find 

(23) 

By composing the evolution in (i) and (ii) we find 



Po(2tt) = ap (ir) 

Pi (2tt) = Pi (tt) + (1 - a)%(n)e ix 



A = 



a 



a(l — a) 



(1 - a)e lx a + (1 - a) 2 e ix 



(24) 



The generating function is obtained by diagonalization 
of Eq. 

2 

A M = a+(l-a)| V(l 1 «)V + 4a+(l- a ) 2 ^- , (25) 

where we introduced the short hand notation y = e %x ^ 2 . 
Current and noise follows by differentiation: 



I 



1 - a 
1 + a 



P = Aa 



1 - a 
(1 + a) 3 



(26) 



For r <c 1 Eq. (|25|l is very accurate and holds for 
< a < 1. It is instructive to study its behavior in the 
two opposite limits of a <C 1 and 1 — a <C 1. 

For a <C 1 the probability of transferring the electron 
during the half cycle is nearly 1 . Linearizing Eq. (|25|l in a 
wc find the generating function of a binomial distribution: 



Sn(x) 



[2a + (1 - 2a)e lx ] 



i X ]N 



(27) 



This means that at each cycle one electron is transmit- 
ted with probability 1 — 2a. The cycles are indepen- 
dent: after N cycles the probability of having trans- 
mitted n electrons is simply given by the binomial dis- 
tribution (1 — 2a) n (2a) N ~ n . Cycles are indepen- 
dent for a — > since at each cycle the system is reset 
to the stationary state within accuracy a 2 , regardless 
of the initial state. The stationary solution is given by 
the eigenvector of Eq. (|24|) with eigenvalue 1 for \ = 0: 
\p.st) = {^/(l + a), 1/(1 + a)}. Calculating the transmis- 
sion probability for one electron during one cycle with 
initial condition given by \p s t) one obtains, with linear 
accuracy, the correct result 1 — 2a appearing in Eq. (|27|l . 

For a — ► 1 the probability for one electron to tunnel 
during a cycle is very small. We can thus expand the 
generating function in the positive quantity (1 — a) <C 1. 
This gives the following surprising result: 



-S N (x) = , 



+ (1 - a)e 



»x/2 



(28) 



We find again that the periodicity of the generating func- 
tion has changed. Eq. (|28|l describles a system of 1/2 
charges that at each cycle have a probability 1 — a of 
being transmitted. The situation is similar to the static 
case. We can again create a mapping on a fictitious sys- 
tem of charges 1/2 and say that every time that one elec- 
tron succeeds in jumping on or off the central island, one 



charge 1/2 is transmitted in the fictitious system. This is 
possible, since it is extremely unlikely that one electron 
can perform the full shuttling in one cycle. Thus after 
many cycles (N ^> 1) the counting statistics of these 
two systems coincide. The cycles are no more indepen- 
dent like in the case for a«l, but the problem can be 
mapped onto an independent tunnelling one. For a inter- 
mediate it is more difficult to give a simple interpretation 
of Eq. Q25[) . since different cycles are correlated in a non 
trivial way. 

In Ref. [l5| the current and noise within a similar model 
have been calculated, but only for a single cycle using 
the the stationary solution as initial condition. This ap- 
proach clearly neglects correlations among different cy- 
cles. We have seen that this is an excellent approxi- 
mation for a — ► 0: Eq. (|27|l represents N uncorrelated 
cycles. But it fails completely in the opposite limit of 
a — ► 1, where the main contribution to the current fluc- 
tuations comes from the cycle-cycle correlations. This 
can be seen as follows. Starting from the stationary so- 
lution for a — > 1 (i.e. {1/2,1/2}) one can calculate the 
average number of particles transmitted n = (1 — a)/2 
and its fluctuation (n — n) 2 = (1 — a)/2 during a single 
cycle. From Eq. I|28|) we see that the average current over 
a large number of cycles is correctly reproduced, but the 
noise differs by a factor of 2. This difference increases 
with higher moments. Even if the fluctuation during a 
single cycle is an interesting physical quantity, the exper- 
imentally relevant one is the long time fluctuations. 



"Mixing" regions 

We expect that Eq. (|25|l describes pretty well the 
counting statistics P/v(") for n < N, but it is clear 
that it fails completely for n > N for which it gives 
Pjv( n > N) = identically. As a matter of fact, the ap- 
proximation does not take into account than more than 
one electron per cycle can be transmitted. This is an 
artefact of the assumption that Tr, and Tr are never 
non-vanishing at the same time. 

In order to improve the approximation, but keeping the 
problem solvable analytically, we need to treat differently 
the left and right "mixing" regions cf> f» 0, 7T, 2tt. In these 
regions at lowest order Tl ~ Tjj. We thus divide the 
time evolutions in 5 steps. For < <^>0i \4> ~ tt] < <f>0i and 
2ir — <f> < <fio ~ 1/a <C 1 wc calculate the evolution with 
T L = Tr = r°. For 0o < 4> < Tr - 0o and n + O < 4> < 
2tt — 0o we use instead the previous approximation for 
regions (i) and (ii). The approximation is summarized in 
Fig. |3 where the exact and the approximate dependence 
of (4>) is shown. 

The contribution to A of the three mixing regions de- 
pends on 0o and T only through their product T0o = r. 
The approximation is meaningful only for a small, oth- 
erwise the constant approximation for the probabilities 
in the mixing regions would not be accurate. We thus 
consider only the small a limit. Keeping linear terms in 
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FIG. 3: Exact and approximate dependence on <j> of Tl used 
to obtain analytically the FCS in the large oscillation limit. 
Tn((j>) has the same form shifted by ir in the <j> axis. 

a we have 

— 4t 

AmW = — [(l-2a)(y 2 - 1) + 2j/ sinh(4rj/) 

+ (1 + 2a + (1 - 2a)y 2 ) cosh(4r?/)] (29) 

with the same short hand notation y = e lx ^ 2 . 

This expression, through —Sn(x)/N = lnAjvf(x) gives 
the FCS for large a in different limits. Let us begin with 
the case r <C 1. Expanding Eq. (|2"§)) up to second order 
in t we obtain 

e -s N (x) = Xm{x) n = ^ Q+/3ie ix + ^ e 2ixf (30) 
with 

ft = 2a(l-4r + 8r 2 ) 

ft = l-2a(l-4r + 4r 2 )-4r 2 (31) 
P2 = 4r 2 (l-2a) . 

The interpretation is again simple, Eq. i|3U|) gives a tri- 
nomial distribution, at each cycle there is a probability 
ft of transmitting ft electrons per cycle. This approx- 
imation holds for any < ft < 2 and t « T/a <C 1, 
a> 1. In this case the probability of transmitting more 
than two electrons within a cycle is extremely small since 
Ph>i ~ 7" 2n— 2 . The importance of the parameter r for 
0n>i proves that to understand the probability of charge 
transfer for n > N it is crucial to correctly treat the 
mixing regions where Tl and Tr are both non- vanishing. 
This means that in this limit the FCS for ft = n/N > 1 is 
determined mainly by the value of r, while for < ft < 1 
is a that controls the FCS. 

Using the saddle point approximation, it is possible to 
obtain explicitly the probability for ft ~ 1: 

ln[P N (n)] = f h A - (1 - ») ^ for ft < 1 

AT " | ]n ft - (ft - 1) In for ft > 1 . 

(32) 

The probability has a sharp maximum at ft = 1, as ex- 
pected, and its logarithm decreases approximately lin- 
early on both sides, with slopes controlled by two differ- 
ent parameters. For ft < 1 the slope is approximately 



given by — ln(2a), while for ft > 1 it is given by ln(4r 2 ). 
Since the parameter a decreases exponentially with a, 
while r is only inversely proportional to a, the peak 
around ft = 1 is asymmetric with an excess to the left 
for moderately large a, and with an excess to the right 
for larger a. 

When r > 1 the previous expansion in r cannot be 
used, but in the small region ft « 1 we can find analyti- 
cally the FCS expanding Xm(x) i n powers of y. In fact, 
it turns out that the saddle point equation l|12fl is solved 
by X = i x with x real and large, thus with y <C 1. We 
find that the expansion of Am contains only even powers 
of y and at the fourth order coincide with Eq. I|3U|) , but 
with the (3 coefficients given by (for large r) 

/3 = 2aer iT 

ft = 4(1 + 2a)e~ 4T T 2 (33) 
ft = f(l + 2a)e- 4 -r 4 . 

For small y it exists a region ft/ft < y 2 < where 
the Piy 2 term dominates the other two terms. We thus 
find again the same behavior of Eq. I|32|l for the probabil- 
ity, but with the coefficients given by Eq. I|33|l . Note that 
now ft > ft, this means that the probability of trans- 
mitting more than one electron per cycle is always larger 
than the probability of transmitting less than one per cy- 
cle. Actually for large a the asymmetry is extreme, the 
slope for ft < 1 is much larger than the slope for ft > 1 
which is moderately positive. 

For ft 3> 1 we cannot expand anymore for small y. 
Since a is not crucial to understand this region we can 
set a = into Eq. 

Am(x) = e~ 4r [y cosh(2rj/) + sinh(2ry))] 2 . (34) 

For large real y we thus find that the generating function 
is that of a static grain active [cfr. Eq. (|15fl ] for a fraction 
4</>o/(27r) of the time: 

-S( X )/N = 40 o r (V*/ 2 - l) . (35) 

IV. GENERAL RESULTS 

The results discussed above can be now compared with 
the numerical results valid for arbitrary values of the am- 
plitude a. These are obtained by solving numerically the 
system of differential equations i|13fl to calculate A. The 
matrix is then diagonalized and the maximum eigenvalue 
in modulus selected. Current and noise are obtained by 
numerical differentiation, while the FCS is obtained by 
solving numerically Eq. (|12fl . 

We begin by discussing current and noise. Fig. 0]shows 
the average number of electrons and its fluctuation for 
different values of T as a function of the amplitude a. We 
first notice the qualitative difference between T smaller 
or larger than 1. In the first case the oscillation of the 
central grain largely increases the current, while in the 
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FIG. 4: Current (plain line) and noise (dashed line) as a func- 
tion of the oscillation amplitude for different values of F. From 
the top-left pannel T — 0.01, 0.1, 1, 5. The dashed light lines 
are obtained with the analytical approximation l|29[l . The 
short lines at a — are the static results given by 1151 , 



second case it reduces it. In both cases for large a the 
current saturates towards one electron per cycle. From 
our previous analysis we know that for large T the sat- 
uration happens only for very large a, when r = 2T/a 
becomes small enough to reduce the contribution of the 
central region. (The choice of a factor 2 into the defini- 
tion of r is arbitrary and it simply improves the accuracy 
of the analytical approximation. Any factor of the order 
of one does not change significantly the results.) A strik- 
ing feature that appears from the plot for the noise is the 
enormous reduction of the Fano factor. The transport 
becomes deterministic due to the shuttling, it is very 
difficult that the grains perform an oscillation without 
transmitting one electron. 

We believe that measuring noise and current in a de- 
vice can give a clean indication if the system is actually 
shuttling electrons. It can discriminate between a sim- 
ple coupling between the mechanical and the electronic 
degrees of freedom of the system not associated with the 
shuttling mechanism. 

In Fig. 0] we also plot the comparison with our sim- 
ple analytical approximation for large and small a. The 
agreement is pretty good, indicating that the crucial fea- 
tures are correctly reproduced by our simple picture of 
evolution in five steps. 

Let us now discuss the counting statistics. In Fig. Owe 
show the evolution of ln[Pjv(n)]/iV when a is increased 
from to 5 by steps of one unit. We show the case V = 0. 1 
that is a good representative of the small T limit. The 
full evolution from the static (a = and I = ttT) to 
the deep shuttling regime (a = 5, I re 1) is obtained 
(cfr. also the current in Fig. 0J. The maximum of the 
distribution moves from approximately 7r/10 to 1. Two 
features are particularly striking: (i) the peak becomes 
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FIG. 5: ln[Pjv(n)]/iV as a function of the number of electrons 
transferred per cycle h for F = 0.1, and oscillation amplitude 
a — 0, 1, 2, 3, 4, and 5 (from left to right). 
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FIG. 6: The same as Fig. EJfor V = 1 



very sharp at the point that a discontinuity of the slope 
of InP appears at n = 1 (it) it becomes asymmetric. The 
fact that the peak is symmetric in the static case is not 
surprising, the probability of transferring more or less 
electrons than the average should not be very different. 
When a becomes large we have instead shown that those 
probabilities are controlled by two different parameters, 
for h < 1 by a and for n > 1 by r = 2T/a. The numerical 
results confirms this prediction. The behavior around the 
maximum is well described by the (nearly linear) form 

Figure [S] shows the case r = 1. In contrast with the 
previous case now for a = the maximum of the dis- 
tribution if for n = 27r, larger than 1. Shuttling will 
reduce the current to 1. The main contribution to the 
transport comes from the sequential hopping through the 
grain when both Tl and Tr are non vanishing. The oscil- 
lation reduces this region in favor of regions where only 
one r is non vanishing. In this limit one electron per 
cycle is transferred. Since this regime is attained when 
the contribution of the region x ~ becomes negligible, 
i.e. when r = 2r/a — > 0, this means that one needs 
huge oscillation amplitudes to reach the truly shuttling 
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regime of h = 1. For large a, but not yet in this limit, 
the probability has the form shown in Fig. We con- 
sidered also this limit analytically after Eq. (J2HJ. Like 
in the previous cases a singularity develops at 1, but in 
this case the probability remains monotonic at 1 (for not 
too large a). The effect of the shuttling is thus mainly 
to enormously reduce the probability that less than one 
electron is transferred, and then to slightly shift the max- 
imum in the distribution from 7rr > 1 towards 1. This 
is due to the fact that due to the oscillations at least 
one particle is always transferred and the probability of 
transferring more than one particle is reduced, since the 
time spent by the shuttle in the central region [cc| <gC A is 
shorter. 



V. CONCLUSIONS 

In conclusion we have studied the full counting statis- 
tics of charge transfer in a single electron transition struc- 
ture where the central grain can oscillate at a given fre- 
quency. The two relevant parameters are the oscillation 
amplitude divided by the scale of the exponential depen- 
dence of the resistance (a), and the probability of a tun- 
nelling event during the time 1/lo for the static structure 
(r). We have obtained both numerical and analytical ex- 
pressions for the FCS. The results apply to both driven 
or self oscillating shuttles, when the fluctuation of the 
amplitude of oscillation can be neglected. The probabil- 
ity of transferring h electrons changes qualitatively as a 
function of a and T. When T > 1 the tunnelling events 



happening when the shuttles passes through the region 
x k, arc always important and very large shuttling am- 
plitudes are necessary to have a well defined shuttling 
regime. 

We also discussed in some details the first two moments 
of the FCS: the current and the noise. We found quanti- 
tative prediction for the reduction of the Fano factor for 
large oscillation amplitudes. 

The study of the FCS permits to understand more 
deeply the dynamics of charge transfer. In some cases 
we found that the effective elementary charge becomes 
1/2 the actual one, due to correlations. This both in 
the static and in the dynamic regime. In other limiting 
cases the statistics is in general polynomial, taking into 
account the probability of different outcomes at each cy- 
cle. Generalization of the theory to a larger number of 
available states in the grain, or the inclusion of an asym- 
metric hopping probability is straightforward and can be 
important to study more realistic systems. 
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